Temperature-pressure scaling for air-fluidized grains on approaches to Point J 
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We present experiments on a monolayer of air-fluidized beads in which a jamming transition 
is approached by increasing pressure, increasing packing fraction, and decreasing kinetic energy. 
This is accomplished, along with a noninvasive measurement of pressure, by tilting the system and 
examining behavior vs depth. We construct an equation of state and analyze relaxation time vs 
effective temperature. By making time and effective temperature dimensionless using factors of 
pressure, bead size, and bead mass, we obtain a good collapse of the data but to a functional form 
that differs from that of thermal hard-sphere systems. The relaxation time appears to diverge only 
as the effective temperature to pressure ratio goes to zero. 



PACS numbers: 64.70.ps, 64.30.-t 

The relaxation time for amorphous liquids can grow 
unbearably long when the temperature is lowered [TJ [5] . 
It can also grow when the pressure is increased, although 
this is more difficult to study experimentally [3]-[5]. Sim- 
ilarly, the relaxation time for colloidal suspensions can 
exceed experimentalists' patience when the packing frac- 
tion or pressure is increased [HHH]- In both the thermal 
and colloidal glass transitions, the particles appear to 
develop a fixed set of neighbors and the bulk medium 
appears to become mechanically rigid. It was recently 
suggested that these two glass transitions are manifesta- 
tions of the same phenomenon for the system of thermal 
hard spheres [S]. In such a system, dimensional analysis 
suggests that the relaxation time r, made dimensionless 
as T(PfT'*~^/m)^/^ by pressure P, the sphere diameter cr, 
and the sphere mass to, must depend only on the dimen- 
sionless ratio T/Pa'^, where d is the dimensionality and 
the Boltzmann constant is set to unity. Thus, the dimen- 
sionless relaxation time increases in exactly the same way 
whether T is lowered or P is raised. 

Although no system behaves exactly like hard spheres, 
Medina-Noyola and coworkers showed that there should 
be "dynamical equivalence," so that soft spheres behave 
as hard spheres with a smaller diameter (TUl HT] . Indeed, 
it was found that particles with a variety of finite-ranged 
repulsive interaction potentials exhibited collapse of di- 
mensionless relaxation time with T /pa'^ as long as the 
pressure was low, so that P ^ e/c'*, where e is the in- 
teraction energy scale [5] . These results suggest that real 
systems might exhibit the temperature-pressure scaling 
expected for hard spheres; however, this has not been 
tested by experiment. For hard-sphere colloids, this is 
not possible for a single system because the packing frac- 
tion (f>, or equivalently the pressure P, is the only con- 
trol parameter; temperature is bounded by the freezing 
and boiling points of the solvent and therefore cannot be 
varied appreciably. For molecular liquids, scaling is not 
expected to hold where van der Waals attractions are 
appreciable compared to the pressure; in this regime the 
scaling must be modified to account for the mean-field 



effect of the attractions [8]. 

Here we describe experiments on a granular monolayer 
of bidisperse beads, subjected to random in-plane forc- 
ing by turbulence in a uniform up-flow of air |12H15| . In 
this system, jamming may be approached by variation 
of either air fiow rate or bead density [12]. The thermal 
energy of these beads is negligible compared with their ki- 
netic energy due to airflow. However, it has been shown 
that there is a well-defined effective temperature, T^s, 
that corresponds equivalently to the kinetic energy, the 
ratio of diffusivity to mobility, and the energy of an em- 
bedded harmonic oscillator [15j. Here, we push the con- 
cept of effective temperature further by asking whether 
temperature- pressure scaling is obeyed using Toff. 

Indeed, we find reasonable collapse with T^g/Pa^ of 
both the equation of state and relaxation time, even 
though our system is far from equilibrium. However, the 
functional dependence on the temperature-pressure ra- 
tio is different from that in a true thermal system. These 
results suggest that while the effective thermal glass tran- 
sition and the colloidal glass transition remain the same 
phenomenon even for these driven hard spheres, the glass 
transition dynamics differ from those of equilibrium hard 
spheres. 

Fluidization and imaging procedures follow Refs. |12| - 
I16| . The granular medium is a 1:1 mixture of 796 steel 
bearings, with diameters a = 0.397 cm and lAa, which 
uniformly fill the (14.9 cm)^ sample square to 67% pro- 
jected area fraction when the system is level. The super- 
ficial airflow speed is 700 cm/s, which is sub-levitating 
and corresponds to Re=2000. To measure pressure, 
we tilt the entire apparatus to induce a component of 
gravity along the monolayer and determine the mass 
per unit length that is above a given depth. When 
tilted by angle 9, ranging here from 0.18° to 0.90°, the 
grains are fluidized near the top but become progres- 
sively jammed with increasing depth z below the top 
edge. Video data are collected at 120 frames per sec- 
ond, typically for 20 minutes, and are analyzed in strips 
of width Az = 0.84 cm. At each depth we determine 
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FIG. 1: (Color online) (tphivp.pdf) Projected-area packing 
fraction, effective temperature, and relaxation time, all plot- 
ted parametrically versus pressure for different tilt angles as 
labeled. 



the time average of the packing fraction 0, the pressure 
P = rrizgsmO / L where is the mass of all beads be- 
tween z and the top edge, and the mean-squared dis- 
placement. From the latter we deduce both the granular 
effective temperature T^s as the average bead kinetic en- 
ergy, and the relaxation time r. Here, t is defined as the 
time needed for the root-mean-squared displacement to 
equal a. Technical details are available on-line |17) . 

The final results for </>, Tog, and r are plotted vs pres- 
sure P in Fig. [T] for five different tilt angles 9. Each point 
corresponds to a different depth, such that the pressure 
is zero at the top edge and increases with depth. Fig. [T] 
shows that as P increases towards the bottom, jamming 
is approached in that both cf) and r increase while Teff 
decreases. Note that at a given P there is a substantial 
range in (/), Tog, and r values for the different tilt angles; 
therefore, any collapse achieved by temperature-pressure 
scaling will have significance. 

The first quantity we consider is the equation of state, 
which is independent of the bead dynamics. Thus we 
form the dimensionless combination Pa^/Tes and plot 
it vs (/) in Fig. [5] This causes the raw data shown in 
Fig. [pL-b to collapse reasonably well onto a single curve, 
with scatter that is random in tilt angle. As expected 
the dimensionless pressure is low for small 0, grows with 
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FIG. 2: (Color online) (eqofstateall.pdf) Equation of state for 
pressure divided by temperature, vs packing fraction, where 
a is the small bead diameter. The solid curve is for a two- 
dimensional simulation of disks interacting via hard-core re- 
pulsion. The short and long dashed curves are fits to free- 
volume and Carnahan-Starling forms, respectively. 



increasing </>, and appears to diverge as (/)—>■ 0c where 
(j)c ~ 0.84 corresponds to random close packing for hard 
spheres in two dimensions [18]. For comparison, we show 
fits to two forms: free-volume theory, Pa^ /T^s oc </>/[! — 
{(j)/(j)cy/^], and Carnahan-Starling [H HD]- Both give 
a reasonable description, but the former fits better near 
jamming and the latter fits better away from jamming. 

We note that the collapse of Pa^ /T^s with is remark- 
able, given the nonequilibrium nature of the system. Our 
result shows that the effective temperature gives rise to 
a well-defined equation of state, underscoring the con- 
clusion of Ref. that the air-fiuidized beads have an 
effective temperature with thermodynamic meaning. 

The measured equation of state may also be compared 
with simulation results for a thermal system. For this, 
event-driven molecular dynamics are performed in the 
microcanonical ensemble for a 1:1 bidisperse mixture of 
1024 hard-core disks with the same diameter and mass 
ratios as in the experiment. The results for Pct^/T are 
plotted as a solid curve vs 4> in Fig. |2] In comparison 
with experimental data, the simulated equation of state 
is slightly high for (j) < 0.5 and does not diverge rapidly 
enough for (p > 0.8. One possible source of this dis- 
crepancy is that the air-mediated bead-bead repulsion is 
strong or long-ranged, as measured previously for two 
beads alone [16]. However, the interactions are based 
on turbulent wakes and hence are not pair-wise additive; 
plus, a similar bead-boundary interaction is washed out 
for many-bead systems. Furthermore, if air-mediated in- 
teractions (pairwise-additive or not) were significant, we 
would not expect to find collapse, which should only oc- 
cur in the hard-sphere limit. Thus we speculate that the 



3 




0.4 0.5 0.6 0.7 0.8 0.9 f ' 

Packing Fraction, (j) 0.01 0.1 J / Po^ 1 



FIG. 3: (Color online) (trelaxvphiA.pdf) Scaled relaxation 
time versus packing fraction, where P is pressure and m is 
the small bead mass, for different tilt angles as labeled. The 
solid curve is for a two-dimensional simulation of disks inter- 
acting via hard-core repulsion. 



nonequilbrium nature of the experimental system is re- 
sponsible for the discrepancy with the simulation. The 
driven system has a well-defined equation of state that 
just happens to differ slightly from the equilibrium equa- 
tion of state. 

Next we consider the relaxation time, t, which is the 
key dynamical quantity specifying the extent to which 
the system is jammed. Since the Reynolds number is 
large and the dynamics are coUisional, the viscosity of air 
is not relevant for setting the time scale. Thus we form 
the dimensionless combination r(P/m)^/^, as in Ref. [S], 
and plot it vs packing fraction in Fig. [3] Data for the dif- 
ferent tilt angles all collapse reasonably well onto a single 
curve. Again, the deviation from collapse is not system- 
atic in tilt angle. Note that t(P/to)^^^ decreases towards 
a number of order unity at low (j) far from jamming, re- 
inforcing the conclusion that T{P/m)^^'^ is indeed the 
correct dimensionless relaxation time to consider. Above 
(j) ~ 0.7, the relaxation time grows at an ever increasing 
rate and appears to diverges as — >■ </)c as expected. The 
simulation results for equilibrium hard spheres, in com- 
parison, are lower at small volume fractions and increase 
more steeply at high volume fractions. 

We also plot the scaled relaxation time rfP/m)^/^ ver- 
sus scaled temperature Teg/Pcr^ in Fig. Wl There we 
also include data from earlier experiments (l2'-'15' , where 
Teff/Pcr^ values were obtained by interpolating the mea- 
sured equation of state in Fig. [2] to the desired packing 
fractions. The good agreement between prior and cur- 
rent data shows that analyzing the video data in nar- 
row strips does not introduce unwanted artifacts. It also 
shows that the bead sizes and masses, which are all dif- 
ferent except for Refs. [El [14], do not noticeably affect 
the collapse. The common behavior of all data, for over 
two decades in scaled temperature and over five decades 



FIG. 4: (Color online) (tauvT.pdf) Scaling plot of relaxation 
time vs temperature, where P is pressure and m and a are 
respectively the mass and diameter of the small beads. The 
solid circles are the same data shown in Fig. 3] The open 
circles are earlier data, compiled from Refs. |12fll5 l and with 
Toff / Pcj^ taken according to Fig. [2] and the known packing 
fractions. The solid curve is for a two-dimensional simulation 
of disks interacting via hard-core repulsion, while the other 
three curves are fits given by Eq. |T| as labeled. 



in relaxation time, is quite different from the simulation 
results for equilibrium hard spheres. In particular, the 
experimental relaxation times appear to diverge only in 
the limit Tcs/Pcr^ — 0, while the simulated relaxation 
times appear to diverge more rapidly. 

In the on-line supplement [I^ we also show that good 
collapse is found for the size of dynamical heterogeneities 
versus the dimensionless relaxation time [211 [22], both for 
air-fluidized beads [TBHT5] as well as for two other systems 
on approach to jamming [251 124| . 

The Fig. [4] scaling plot of relaxation time versus tem- 
perature allows comparison with well-known functions 
for thermal systems. For this we display fits to the 
Vogel-Fulcher-Tammann (VFT) and Elmatad-Chandler- 
Garrahan (EGG) [25] forms, as well as a stretched expo- 
nential, respectively 



'(P/to) 



r exp[A/(a; - x^)], 
cx I exp[B(l/a; - l/xif] 
Uxp[C/a;"], 



(1) 



where x = Tcs/Pa^ and the other quantities are fitting 
parameters. The VFT form traditionally accounts for 
divergence of relaxation time at a nonzero temperature. 
But here the best fit gives a proportionality constant of 
1.3 ± 0.2, A — 0.36 ± 0.02, and a negative critical tem- 
perature Xo — —0.029 ± 0.005. The latter causes the 
fit to roll over to a constant at T^s/Pa'^ goes to zero, 
which is not a physical feature supported by the data. 
The best fit to the EGG form gives a proportionality 
constant of 0.0035 ± 0.0025, B = 0.00032 ± 0.00002, and 
= -0.00625 ± 0.00015. This diverges at zero tem- 
perature, and gives a better fit than VFT; however, the 
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divergence is too fast. Overall, the best fit is a stretched 
exponential with proportionality constant 0.40 ± 0.05, 
C = 2.0 ± 0.6, and stretching exponent a — 0.35 ± 0.04. 
Thus it appears that the data diverge only at zero tem- 
perature. 

In the case where T^e is well-defined so that different 
definitions yield the same result, one would expect col- 
lapse with Tos / Po^ with the same functional form as in 
equilibrium. Indeed, simulations of thermal hard spheres 
under shear |26j support this expectation. It is surpris- 
ing that the relaxation time for our driven system shows 
temperature-pressure scaling but with a functional form 
that is different from that in equilibrium. The scaling 
collapse of relaxation time in our system suggests that 
the effective thermal glass transition, observed by low- 
ering Tcff, is equivalent to the "colloidal" glass transi- 
tion, observed by varying P. However, the difference in 
functional form suggests that the glass transitions of the 
driven system are somewhat different from those of the 
equilibrium system. In particular, for equilibrium hard 



spheres, the value of T jPa'^ at which the relaxation time 
diverges is ambiguous; the VFT form, which diverges at 
a nonzero value of TlPa'^, and the ECG form, which di- 
verges at TjPa'^ = 0, both fit simulation data equally 
well It is still not known whether equilibrium hard 
spheres have a thermodynamic glass transition (corre- 
sponding to a divergence at nonzero T/Pa'^). But for 
the driven granular hard spheres, the answer seems much 
more clear: The relaxation time appears to diverge at 
T/Pa'^ = 0. This has a special significance because it 
corresponds to the zero-temperature jamming transition 
of spheres. Point J [37|[2S]. Our results therefore imply 
that glassy dynamics in air-fluidized grains are controlled 
not by a thermodynamic glass transition, but by Point J. 
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